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^ : ABSTRACT 
pu 

According to the core accretion theory, circumbinary embryos can form only 
beyond a critical semimajor axis (CSMA). However, due to the relatively high 
density of solid materials in the inner disk, significant amount of small planetesi- 
mals must exist in the inner zone when embryos were forming outside this CSMA. 
So embryos migration induced by the planetesimal swarm is possible after the 
gas disk depletion. Through numerical simulations, we found (i) the scattering- 
driven inward migration of embryos is robust, planets can form in the habitable 
zone if we adopt a mass distribution of MMSN-like disk; (ii) the total mass of 
the planetesimals in the inner region and continuous embryo-embryo scattering 
are two key factors that cause significant embryo migrations; (iii) the scattering- 
driven migration of embryos is a natural water-deliver mechanism. We propose 
that planet detections should focus on the close binary with its habitable zone 
near CSMA. 

Subject headings: planetary systems — binaries: close — planets and satellites: 
formation 



1. Introduction 



Circumbinary planets - planets encircle both members of a binary are interesting celes- 
tial bodies discovered by exoplanet-hunting missions. So far, six exopla nets circling main - 
sequence binary stars were found by Kepler mission, they are Kepler-16 b (IDoyle et al.ll201lh . 
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Kepler-34 b, Ke Dler-35 b flWelsh et al.ll2012h . Kepler-38 b flOrosz et al.ll2012al ). Kepler-47 b 
and Kepler-47 c (jOrosz et al.ll2012bl ). Though exoplanet-hunting plans selectively avoid close 
binaries, circumbinary planets may be common. Researches manifest that more than 1% of 
close binary stars have giant planets in nearly coplanar orbits, which yielding a Galactic 
population of at least several million ( IWelsh et al.l 120121 ) . All the six circumbinary planets 
discovered hitherto by Kepler are Saturn-like or Neptune-like planets, they are not suitable 
for the existence of extraterrestrial life. We know the ultimate aim of Kepler is to find 
Earth-like exoplanets, and to search for extrasolar life. It is natural to ask whether terres- 
trial planets can form in the circumbinary case? This Letter intends to explore some aspects 
of this topic. 

The semimajor axis (SMA), eccentricity and mass ratio of a binary are designated by 
Ob , cb and n = MbI{Ma+M b) {Ma is the bigger one) in the paper, respectively. According 
to iHolman fc Wiegertl (119991 ). there is a critical SMA (od) within which test particles are 



unstable (beside some specific initial configurations). Based on the core accretion model, 
there is another critical SMA (ac2, usually outside of Od) within which planetesimal accretion 
is impossible due to the strong perturbation of the binary. In this Letter, we propose a 
mechanism to form terrestrial planets in the habitable zone (HZ) through inward migration 
of embryos outside ac2- It is organized as follows: a summary of former works and our ideas 
are presented in Section 2; Session 3 and 4 show our model and simulation results. The 
conclusions are summarized in Section 5. 



2. Synopsis of former works and our ideas 



From planetsimals to embryos. In the circumbinary case, the perturbations of the com- 
panion star can stimulate the eccentricities of the planetesim als, which lead to high i mpact 
velocities between planetesimals. By ignoring the effect of gas, iMoriwaki fc Nakagawal (120041 ) 
concluded that the planetesimal accretion cannot occur within 13 AU {n = 0.2, = 1, 
cb = 0.1). The coupled effects of companion perturbations and gas friction lead to strong 
size- dependent orbital pha sing, which reduce s the planetesimal impact velocities and makes 
accretion more optimistic (jScholl et al.ll2007l ). However, when the disk evolution is included, 
pericenters of planetesimals are less collimated than that in t he stationary and ax isymmetric 
disk model, which again inhibits the planetesimal accretion (IMarzari et al.ll2008l ). 



As a result, iMeschiaril (j2012aj) con cluded pla r ietesira al accretion might be effectively 
inhibited within 4 AU in Kepler-1 6 (AB). IMeschiaril (l2012bl ) considered the turbulence driven 
by the magnetorotation instability. He fo und that planetesinial acc retion can be inhibited 
even in 4-10 AU region in Kepler-16(AB). iPaardekooper et al.l (120121 ) confirmed that in situ 
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growth starting from planetesimals smaller than 
34b, and Kepler 35b. 



10 km is difficult for Kepler 16b, Kepler 



From the embryos to the terrestrial planets. I Quint ana fc Lissauerl (j2006[ ) found the 
apastron distance Qb = clb{^ + ^b) is a critical parameter for the formation of circumbinary 
terrestrial planets. When Qb < 0.2 AU (for 0.05 AU < as < 0.4 AU), the formed planetary 
systems are very similar to those around single stars, whereas those with larger Qb tend to 
harbour fewer planets, especially interior to 1 AU. However, they didn't consider the origin 
of the initial embryos. 

Since circumb inary embryos ca n only form in large distances (beyond the ice line in 
Kepler-16 system (IMeschiaril l2012al )). embryos' inward migration is the only path to form 
terrestrial planets i n the inner zone. There are two kinds o f migrations: type I migration 



driven by the planetesimal disk (llda et al. 



driven by gas disk ( IGoldreich fc Tremaine Il979l : IWardI 119971) and the migration of embryos 



2000h . IPierens fc Nelson! J2007[ boosi ) have shown 



that planets with terrestrial to Saturn mass always migrate inwards till they reach the inner 
border of the disk truncated by the tidal force of the binary. However, from the viewpoint of 
planet formation, there are some inevitable problems for the formation of the circumbinary 
embryos or protoplanets. 



Long accretion timescales of distant embryos. Outside ac2, planetesimal accretion will 
be slower than that in single-star systems due to the perturbations of the companion 
star, with the so-called type II runaway growth faortenkamp et ah 2001). On the other 
hand, is a function of time, over longer timescales, the impa ct velocities in t he outer 
regions will also increase, which makes accretion more slowly ( iThebaultl l2006l ) . 



• Rap id dissipation of the circ umbinary gas disk. Simulations of lArtymowiczl (119941 ) 
and iFragner fc Nelson I (120101 ) manifest that the disk will be truncated, distorted and 
teared quickly, so the the lifetime of the circumbinary gas disk will be shorter. Some 
evidences for the reduced lifetimes (~ 0.1-1 Myr) of t he primordial circu mstellar disk 
in binary system had been conffimed by observations (jCieza et al.ll2009l ). 



• Uncertainty of the type I migration. Recent studie s show that the orien t ation and m i- 
gration rate of type I migration is very uncertain (IKley fc Nelson! 120121 : lOrmell |2012| ). 
Inward migration can not be guaranteed in all kinds of disk environments. A notable 
exarnple is the solar system where planets might not have obvious gas-driven migra- 
tion!^ On the other hand, type I migration is unimportant for the terrestrial planets 



new point can be found in IWalsh et al.l ( 20121 ). 
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(iBromley &: KenyonI l201ll ). The reason is that building blocks of solid planets are 
packed too closely to migrate. 



Long accretion timescales of embryos and the reduced lifetime of gas disks indicate that 
significant gas-disk migration of circumbinary planet may be less significant as in single star 
systems. However, due to the relatively high density of solid materials in the inner disk, 
significant amount of small planetesimals must exist in the inner region when small embryos 
was forming outside of ac2- Thus planetesimal-driven migration of embryos maybe work. In 
this paper, we consider the migration and accretion of embryos in a planetesimal disk — a 
'gas free' case that the influence of the residual gas disk can be ignore d. Our work can 
be thought as a complementary work of iPierens fc NelsonI (120071 . 120081 ) — a case that gas 
disk has dissipated when the embryos formed outside of ac2- The aim is to check whether 
the terrestrial planets c an form in the inne r region such as habitable zone (HZ) under such 
'gas free' circumstance. iPayne et al.l ( 120091 ) have discussed the migration of S-type embryos 
in a planetesimal disk, we apply a similar scattering-driven migration mechanism in the 
circumbinary case. 



Numerical Model and Initial Conditions 



The parameter space of the binary systems is relatively large. Here, we select the Kepler- 
16(AB) as the fiducial binary configuration. The fitting formula found in lHolman fc Wiegert 
( 119991 ) is used to calculate the unstable boundary a^i (error bar is ignored). 



a^i = [1.60 + S.les - 2.22e| + 4.12/i- 
4.27eByU - 5.09/^2 + 4.61e|/i2] as- 



(1) 



Using the parameter o: 



Keple r-IGCAB) (see Table 1), we get a^i ^ 0.66 AU. The = 4 AU 



derived by iMeschiaril (j2012al ) is used to define the boundary inside which the embryos for- 



mation is effectively inhibited. We call a^i to the inner regior i and > a^? the outer region. 
According to the 'minimum-mass solar nebula' (MMSN) disk (jHayashilll98ll ). the mass of 
solids within the inner region is about 6-8 M® (scaled by (M^ -|- Mb)/Mq). We suppose 
the total mass of planetesimals in the inner region is 6 in our model. Planetesimals are 
randomly distributed according to the a~" profile. 300 planetesimals are put in the inner 
region, each with a mass of 0.02 M®. Outside ac2, we assume the runaway growth has been 
taken place, resulting a chain of isolated embryos. Because the final r nass profile of e r nbryo s 
is unclear in the circumbinary case, we take the a-disk model given by lKokubo fc Ida! (120021 ) . 
and explore six kinds of possible value of a (a = 0.0, 0.5, 1.0, 1.5, 2.0, 2.5). The distance 
between embryos is 7 mutual Hill radius. We put 100 embryos in the outer region following 
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Kokubo fc Idal (|2002|), namely, 



(2) 



th ^ 0.69 X 10 



-2 1 b \^/'(^^)V2 



10 



X 



(lAu) 



(l/2)(4-a)/M. A"^/^ 



AU, 



(3) 



where b is the orbital separation of embryos (scaled by the mutual Hill radius). Si is the 
reference surface density at 1 AU, fice is the ice factor and is the mass of the central 
stars. Because a ll the embryos lie outside the ice line of the Kepler-16(AB) (~2.3 AU) 
(IMeschiaril l2012al ). fice'^i = ^ice can be treated as a single parameter. As we only focus 
on the formation of Earth-like planets, the mass of embryos is limited between 0.01 to 0.1 
(Moon-like to Mars-like mass). It can be obtained by carefully adjusting the parameter 
Sjce, resulting the total mass of embryos ~ 6M0 in the outer region (4.0 AU to 6.6 AU). 
We found the total mass of embryos in this region can be compatible with the MMSN disk 
(~ in this region). 

We integrate the system to 10^ yr. The Bulirsch-Stoer integrator in MERCURY pack- 
age (j Chambers! Il999l ) is modified to model the close binary configuration. The accuracy 
parameter of 10~^^ is used, all integrations conserved total energy and angular momentum 
within 10"^. 



4. Simulation Results and Analysis 

4.1. The Fiducial Case 

We take an MMSN disk {a = 1.5) as the fiducial case. Planetesimals are simulated as 
test particles (with mass), whereas embryos are modeled as fully interacting bodies. Because 
scattering is a chaotic process, we perform 20 simulations for such a disk. The angular orbital 
elements of planetesimals and embryos are randomly selected. A representative outcome of 
the fiducial case is given in Fig. 1. After 1 Myr, the original 3 embryos outside have 
migrated into the inner region. The innermost one has a SMA of 0.95 AU, the mass is 
= 0.89 M^, and the eccentricity is 0.037. To compare with different runs, we define a 
migration rate as 
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where aimj is the final SMA of the innermost embryo and aim,i is the initial SMA of the 
same embryo. We find r] = — O.TSI^q qI for the 20 runs, where the errors denote the upper 
and lower quart iles. 

In Fig. 2, we give migration rate as a function of time, averaged over 20 simulations. 
We find that, after 2 x 10^ yr (r/ = -OmtoZ), 

no significant migration takes place, which 
indicates that the system has reached an equilibrium state, thus the statistical results are 
reliable. In the fiducially model, we find 2 out of 20 runs (so, 10% chance) result in a planet 
in the HZ of Kepler-16(AB) (Table 1). 



4.2. Parameter Exploration 

Total mass of embryos (Mtot)- We vary the number of embryos outside and keeping 
other parameters unchanged. Ai"e=200, 150, 100, 50, 25 are explored. Moreover, an extreme 
case A^e=l is also studied. We conduct 20 runs for each N^, results are shown in Table 2. We 
find Ne has no meaningful effect on the migration rate and the mass of innermost embryo as 
long as > 100. But when A^ < 25, migration rate decreases significantly. At the extreme 
case A'e = 1, embryo has negligible migration. 

Location of CSMA (0^2)- The initial size and distribution of planetesimals, the value 
of the critical accretion velocity, the growth of planetesimals and oth er factors have great 



uncertainty in determinin g a^2- Different mode l s give different ac2- E.g.. lMarzari et al.l (120121 ) 



shows ac2 ~ 3 AU, while IPaardekooper et al.l ( ]2012[ ) find ac2 > 4.4 AU for Kepler 16(AB). 



Here, ac2 is varied to check its effect on migration rate. We consider ac2=[3.0, 3.5, 4.0, 4.5, 
5.0] AU. Sjce is adjusted to keep the total mass of embryos and planetesimals unchanged 
(equal to the fiducial case). The results are shown in Fig. 3a, we find no significant changes 
in the migration rate and the masses of the innermost planets that formed. 

Density profile (a). The material distribution of a circumbinary disk is unclear. We 
consider other five kinds of density profiles: a =0.0, 0.5, 1.0, 2.0, 2.5. The results are shown 
in Fig. 3b. Again, no significant change in the migration rate and mass, which manifests 
that the inward migration of embryos is generally robust. 

Total mass of planetesimals. Here we consider two extreme cases only: 1, no planetesi- 
mals in the inner region; 2, the total mass of planetesimals is double that of fiducial case (ice 
line is in the inner region). Other parameters are unchanged. We find that in the absence 
of planetesimals, migration of the embryos are severely suppressed, generally with Oimj >3 
AU. For example, if we remove the planetesimals in Fig. 1, the SMA of innermost planet is 
~3.5 AU (see right panel of Fig. 1). It shows that significant (inward) migration is caused 
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not only by the scattering between embryos but also by planetesimals disk. If the mass of 
planetesimals is doubled, we have rj = —O.TG^q^I, aim = 0.90 AU. Increasing the total mass 
of the planetesimals causes stronger inward migration (compared to the fiducial case). The 
probability of forming a HZ planet is also increased (20%). 

The number of planetesimals (Nt). In order to improve the computational efficiency, 
we model the planetesimals as test particles with associated mass. To validate that this is a 
reasonable approximation, we perform additional two suits of runs with Nt=500, 1000, the 
total mass of the planetesimas is unchanged. The results are rj = — 0.74^q;q5; t] = — 0.73^q;q4, 
respectively. Obviously, using a larger number of planetesimals gives us the similar results. 



4.3. Other Systems 



We explore the embryos migration in the other two binary configurations: Kepler- 
34(AB) and Kepler-35(AB). We use a sample formula to estimate the boundary of circumbi- 
nary HZ, 



[HZin, HZ 



La + L 



B 



out I 



La 



[0.95, 2.0] AU. 



(5) 







[0.95, 2.0] AU is t he extended HZ of the sun flKasting et all Il993l : iMischna et all I2OOOI : 
Ouarles et al.l 1201 21 ). La + Lb and Lq are luminosi ties of the binary and th e sun. By 



including planetesimal formation and dust accretion, iPaardekooper et al.l (120121 ) find ac2= 
2.76 AU for Kpeler-34(AB) and 0^2=2.7 AU for Kepler-35(AB). The location of t he ice fine 
for an irradiated disk can be estimated from the scaling 2.7[(M^ + Mb) /Mof AU ( Meschiari 
2012ah . We get 10 AU for Kepler-34(AB), and 7.8 AU for Kepler-35(AB), so we can ignore 
the material enhancement caused by the ice line in putting embryos due to Ojce > ac2- 
Other initial conditions have an analogy to the fiducial case of Kepler-16(AB). Results are 
summarized in Table 1. Again, we find that formation of planets in HZ is possible. 



4.4. The effects of a pre-existing giant planet 



Three giant planets discovered in Ke pler-16(AB), -34( AB) and -35(AB) lie just outside 
the stability boundaries of the binaries (IMeschiaril l2012al ). In order to check how a pre- 
existing planet might affect the previous migration scenario, we performed several additional 
simulations by initially placing a planet in orbits similar to Kepler-16 b, -34 b, -35 b. Other 
initial conditions are similar to the fiducial case. We find rj = —0.70~^q'qI and a^m = 1-15 AU 



in Kepler-16(AB). No innermost planets can always stay in HZ (with ap{l + Cp) > HZ 



out I 
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But in Kepler-34(AB), we get rj = -0.44lo;o3 and amin = 1-57 AU, 75% innermost planet 
are located in HZ. Kepler-35(AB) have about 90% innermost planet in HZ. Giant planet 
and binary compose a three-body system, it has a new stable boundary (larger than 
ad), the planetesimal disk that we used in the fiducial case is truncated by this Ocs- As a 
result, the embryo migration is limited by Ocs. Actually, if no giant planets are included in 
Kepler-34(AB) and -35(AB), some innermost planets migrate over HZ, so the fraction of the 
innermost planet in the HZ (Table 1) is smaller than the case with a giant planet. 

5. Summciry 

Inward migration of the circumbinary embryos driven by planetesimal swarm is dis- 
cussed in this Letter, the main results and their implications arc summarized as follows, (i) 
Scattering-driven inward migration of embryos is robust in the circumbinary case. In the 
three close binary configurations we explored, planets can form in the HZ if we adopt mass 
distributions of MMSN-hke disk, (ii) The total mass of the planetesimals in the inner region 
is the main factor that affects the migration (rate). Embryo-embryo scattering is another 
necessary factor that causes significant migration. Migration of a single embryo is limited 
because its increasing mass will reduce the migration rate, ultimately, it clears the adjacent 
planetesimals, migration stops. While in a multiple-embryo system, embryo-embryo scatter- 
ing can continually push the smaller embryo inward, (iii) We propose that planet detection 
should focus on the close binary with its HZ near ac2- (iv) Scattering-driven migration of 
embryos is a nature water-deliver mechanism. In the close binary like Kepler-16(AB), em- 
bryos can only form outside of the snow line (ac2 > Ojce)) so the planet found in the HZ 
should have water contents. 
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Table 1: Binaries parameters and a summary of simulation results. 



Binary 


Kepler- 16 (AB)^ 


Kepler-34(AB)'^ 


Kepler-35(AB)i' 


Ma (Mo) 


0.6897 


1.0479 


0.8877 


Mb (Mq) 


0.2026 


1.0208 


0.8094 


as (AU) 


0.22431 


0.22882 


0.17617 




0.15944 


0.52087 


0.1421 


La + Lb (Lq) 




2.77 


1.35 


aci (AU) 


0.66 


0.84 


0.50 


ae2 (AU) 


4.00= 


2.76^ 


2.70'! 


ttice (AU) 


2.3 


10.0 


7.8 


HZ (AU) 


0.36 ~ 1.02^ 


1.58 ~ 3.33 


1.10 ~ 2.32 




' 3-0.04 


_0 49+0-07 


-0 46+°-^^ 

'-'•^"-0.09 


(Me) 


U.JO_Q 20 


1 iq+0.37 

-'-•-'-^-0.32 


1 1O+0.63 
-'-■■'-O-0.28 


C^mm (AU) 


0.95 


1.33 


1.09 


EIHZ (Fraction) 


10% 


35% 


20% 



" IPovle et al.l |201ll ): ' iWelsh et all |2012h : ' iMeschiaril |2012al) : ^Paardekooper et all |2012l ): ' iOuarles et al 
(I2ni2h 



Note. — amin is the minimum final SMA of the embryos in 20 runs. EIHZ means the fraction of the 
innermost planet in the habitable zone. 



Table 2: Numerical results of changing the number of embryos. 





200 


150 


100 


50 


25 


1 



Mtot 


(Me) 


14.72 


9.81 


5.87 


2.65 


1.27 


0.05 


V 




-0 76+0-15 
"J-o.io 


-0 7fi+0-io 

' 3-0.07 


-0 7Pi+0-05 

"J- ' 3-0.04 


-0 62+0-06 

U-O^-0.03 


-0 tsn+0 07 


-0 06+0-06 

U-UU-0.02 




(Me) 


1 06+1-14 

±.UU_Q 20 


1 ce:+0.99 
l-"3_o.65 


n Q8+115 
U.»0_Q 20 


n qi+0.35 

'-'■^-'--0.07 


1 14+0.11 


4^+0-02 

"-'■^3_o.io 


(^min 


(AU) 


0.93 


0.94 


0.95 


1.39 


1.57 


3.36 
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Planetesimals induced migration 



No planetesimals 



r Myr | 


r 1 0.01 My 


^^Earth scale 


■ 1 ' 


[ I 


* * 


r 0.02 Myr | 


r 1 0.05 My 


o 

• 


- O- 


L 

: 1 f 


o 

# I* 

o • ^ 

• ••: 


t 0.1 Myr 1 ^ 

^ On ° 1 ° 

r o ° ° J° 

o o ° J* 


r 1 0.2 Myr 
■ 1^ * 


o - 

1 • 

•1 

" !* ° • - 

r 


t 0.5 Myr | 


r [ 1 Myr 


• 

• 



Myr 



0.01 Myr 



0.1 Myr 



Myr 



BEarth scale 



3 4 5 
(AU) 



7 1 



3 4 5 
(AU) 



1 2 3 4 5 6 7 
(AU) 



Fig. 1. — The time evolution of the circumbinary planetsimals and embryos. The binary is 
the Kepler-16(AB) configuration (see Table 1). The embryos are represented by filled circles 
whose size are proportional to the physics sizes of the bodies (in scale relative to the Earth). 
The planetesimals are represented by open circles but not in scale relative to the Earth. 
The dashed red line denotes Od and the dashed blue line denotes ac2- The green region is 
the habitable zone (0.36 AU < (1 — e^) and Op (1 + Cp) < 1.02 AU). Right panel: time 
evolution of the circumbinary embryos if the initial planetesimal swarm is removed. 
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Fig. 2. — The value of r] = {aimj — aim, i)/aim, i as a function of time, averaged over 20 runs. 
The median value is plotted, as well as the lower and upper quartile bounds. The inline 
image is zoom of the first 0.2 Myr. 
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Fig. 3. — Variation of the final t] 



"im, I 



)/a. 



im, I 



and the mass of the innermost 



embryos due to changes of the location of (Panel a) or the surface density profile S oc a~° 
(Panel b). Plotted are the median values at 1 Myr (20 runs for each point), the error bars 
are the upper and lower quartiles. 



